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Abstract 

The quantum capacity of thermal noise channel is studied. The extremal input state is obtained at the 
postulation that the coherent information is convex or concave at its vicinity. When the input energy tends 
to infinitive, it is verified by perturbation theory that the coherent information reaches its maximum at the 
product of identical thermal state input. The quantum capacity is obtained for lower noise channel and it 
is equal the one shot capacity. 

One of the most important issues of classical information theory is the Shannon formula, which is the capacity 
of an additive white Gaussian noise channel. It is achieved when the input is a Gaussian noise source with power 
constraint ^ [5]. 

C = l0g2(l + |), (1) 

where S is the power of the source and N is the power of the noise, the bandwidth W should be multiplied 
when it is considered. The formula has guided the design of the practical communication system for decades. 
Correspondingly, in quantum information theory, although a lot of works have been done [2], such a formula 
is remained to be discovered. The Shannon formula comes from the Shannon noisy coding theorem, the later 
gives the capacity of any noisy channel: 

sup I{X;Y), (2) 

X 

where the supremum is taken over all inputs X, I{X;Y) is the Shannon mutual information and Y is the 
output. The counterpart of mutual information in quantum information theory is the coherent information (CI) 
Ic{p,S) = S{£{p)) — S{p^^ ) 0] [Sj. Here S{q) = — Trglog2 Q is the von Neumann entropy, p is the input state, 
the application of the channel £ resulting the output state £{p); p^'^ — {£ (!^T){\tp) {tp\), is the purification 
of the input state p. The quantum channel capacity is 

lim sup -/,(/>„, (3) 

n^oo n 

The righthand side of above formula was firstly proved to be the upper bound of quantum channel capacity p] . 
The equality was proved at the postulation of hashing inequality ^3 ; the hashing inequality was lately established 
[7j. The quantum capacity of a noisy quantum channel is the maximum rate at which coherent information 
can be transmitted through the channel and recovered with arbitrarily good fidelity. For quantum information 
channel can be supplemented by one- or two-way classical channel, thus quantum capacities should be defined 
with these supplementary resources. We here deal with the quantum capacity without any supplementary 
classical channel '3'. 

Quantum capacity exhibits a kind of nonadditivity [H] that makes it extremely hard to deal with. Until 
now, quantum capacity has not been carried out except for quantum erasure channeljHI. We in this paper will 
deal with the quantum capacity of thermal noise quantum channel (which is addressed as Gaussian quantum 
channel before JOI)- The general description of the channel is to map the state p to another state £{p), where 
f is a trace preserving completely positive map. The map £ has a Krauss operator sum representation. That is 
£{p) = AipA\ with ^l^i = I- For additive quantum Gaussian channel, it is quite simple to choose Aa 
to be proportional to the displacement operator ^I] ^U] D {a) = exp[Q;a^ — a*a\. The output state will be 

£{P) - ^ / d^aeM-\a? /N)D{a)pD\a). (4) 
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for the simplest situation of thermal noise channel, where N is the average photon number of the output state 
if the input is the vacuum. 

In dealing with the maximization of the CI, there is a useful lemma in classical information theory which 
gives necessary and sufficient conditions for the global maximum of a convex function of probability distributions 
in terms of the first partial derivatives. The lemma was extended to quantum information theory |12| in 
evaluating the capacities of bosonic Gaussian channels. Let be a convex function on the set of density 
operators which contains /0(o) and p, the necessary and sufficient condition for F achieves maximum on /9(o) is 
that the convex function F{{1 — i)/0(o) + tp) of the real variable t achieves maximum at t = for any p. That is 
^ |t=o-P((l — t)P(o) + ^p) < Generally speaking, CI Ic{p,£) is not a global convex function of its input state 
p. Without lose of generality, let us suppose it is convex lli-ilat the vicinity of some p(o) , the necessary and 
sufficient condition that p(o) is the maximal state will be 



dt 



(5) 



where p is at the vicinity of p(o). The derivative will be— rr(£(p)— £(p(o))) log £{p(^^)) + Tr{p^^ ~ P^o) ) ^^^P^^ ■ 
If £ is a trace preserving completely positive Gaussian operation, then for a gaussian input state p(o) ; the output 
state £(p(o)) ^'^d the joint state p^^ will be Gaussian. Hence their logarithms are quadratic polynomials in 
the corresponding canonical variables|12|. The derivative will be zero under the constraints of the first and 
second moments. Where the trace preserving property of £ is also used. The conclusion is that for input states 
with the same first and second moments, Gaussian input state achieves the maximum of CI for a given trace 
preserving completely positive Gaussian channel. The same conclusion can be applied to CI Ic{pn,£®"), the 
maximum will be arrived when the input state is Gaussian. 

We then turn to Gaussian input state. Every operators A e B{TL) is completely determined by its 
quantum characteristic function xa{z) := rr[AM^(z)], where = exp[— zz-^i?] are Weyl operators and 

R = (Xi, Pi, X2, • • • , P„), with [Xk^ Pi] = iSki- The density operator pn is called Gaussian, if its characteristic 
function xpi^) ^''^^ the form Xp(^) = exp(z?7"^z — ^z^^z). One can show that the first moment i] =Trp„i?,the 
second moment 7 = 2Tr(i? — rj)pn{R — 77)-^ + iJm where 



k=l 



-1 

1 



(6) 



The 2n x 2?! real symmetric matrix 7 is usually called the correlation matrix (CM) of the state p„. The first 
moment represents the displacement of the state, it is irrelative to the problem of entanglement as well as channel 
capacity so that dropped. The completely positive map on the input state pn will be p^'^ = (£(g)I)( |-0„) (ip^l ). 
The CM of the Schmidt purification | ■(/'„) is [H] 



7 f3 
1 



(7) 



where (3 = —0^ = Jn\l —{Jn ~ are purely off-diagonal. It should be noticed that the symplectic 
eigenvalues of 7^ are the square root of the eigenvalues of — (^2x\i7'/j)^' where J2xn — Jn ffi {~Jn) are chosen 
to produce the off-diagonal fi. The application of trace preserving Gaussian channel will result the state p^'^ 
with CM dU 

" Mj7M„ + MlP 
/3^M„ 7 



(8) 



where M„ = M®", = A^f". 

For thermal noise quantum channel, we have Mi = I2, iVi = 2NI2. Then 7^ will be 



7 + 2iVl2„ 13 

1 



(9) 
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The energy of the input state pn is En — J^ii^i + 1/2) =Trp„ J^ii^l^i + 1/2) = (Tr7)/4, where we set the unit 
of the energy such that hu! — 1. Under the energy constraint Tr7 = 4i?„, which state will achieve the maximum 
of the CI? We suppose the state having maximum CI is jO„(o) with its CM 7„(o) = 2£'l2„, where E = En/n is 
the average energy of each mode of input state. We need to prove that 

-Tr((p;)-(p;(o))) logp;(o) + Tr(p^«' - P%)\ogp% < 0, (10) 

where p'^ = f ®"(p„), p^^g^ = ^®"(p„(o))- The density operator p„(o) now is the direct product of pi(o), p„(o) = 
pf(^), we have p'^^^^ = p'^^^^ and p% = pjg®" as well. The CM of p'^^^^ now is 7„(o) + 2iVl2„ = 2(E + 

7V)l2„, thus = {g)^(l - v')!^''' is a thermal state with v' = (E + N - l/2)/(E + N + 1/2). We have 

-Tr((p;)-(p',(o))) logp',(o) = -Tr((p;)-(p;(o))) J2i ^U'' log^'' = under the energy constraint Trp„ + 
1/2) = En, wlicro Ti-p'n E,(4°'J + l/2) = |Tr(7 + 2Afl2„) and Trp^ =Trp'j(o) 1 are used. The density operator 
Pi{p) ^^'^^^ be diagonalized by some unitary transformation [/i, the corresponding symplectic transformation 



cosh rl2 — sinh rJ 
sinh rJ cosh rl2 



(11) 



will diagonalize the CM of pf^.^ , that is <S'i7i(o)5'f = 7i(o)) meanwhile 5i J® (— J)5f = J® (—J). Here tanh2r = 



— 1/(2E + N). The diagonalized CM 71(0) = diag{'jA,jA,7B,'yB}- The density operator p^^>^ is a direct 

product of pf^^ . The unitary transformation diagonalizes p^^-^ will be f7„ = J7®". The corresponding symplectic 
transformation will be 

coshrl2n — sinh rJ„ 
sinh rJn cosh rl2n 



Sn — 



(12) 



Thus UnPn^)Un = Pab' where p^s = (1 — va)v'^°'® (1 — ub)^^'' is a thermal state, with 6, 6^' being the 
annihilation and creation operators of 'reference' R system which is introduced in the purification, and Vj = 
(7i - l)/(7, + 1), U = A, B). Hence Tr(p««' - P%)\ogp% =Tr([/„p«Q'[/t - p%) logp®^ =Tr(C/„p«Q'[/t - 
Pab) log t;^ + log ub)- After the unitary transformation, the density operator UnPn^ is an oper- 

ator function of the creation and annihilation operators ai,ai,b\ and h\. The CM of density operator UnPn^ 
will be Sn'y'^S^, denote it as 

lAA lAB 
IBA IBB 



lU = 



(13) 



with ^AA — cosh^ r(7 + 2iVl2n) — sinh^ rJnjJn + sinh r cosh r(/3J„ — Jn0^), Jbb = cosh^ — sinh^ rJn{-y + 
2Nl2n)Jn + sinhr cosh r(J„/J — 0^ Jn)- Prom the definition of CM, one can get that TrU'„p^'3'f/^ ^.((aja, + 
1/2) \ogVA+{^\hi+l/2) \ogVB) = (Tr7AA log?;A+Tr7BB logt;B)/4, Tt^aa = cosh2r(Tr7)-|-4nA/'cosh^ r - sinh2r 



Tr^-(Jn S)^ - l2n, Tr7Bs = cosh2r(Tr7) + 4niVsinhV - smh2rTryJ-{Jn'^^)^ - hn- When the energy is 
constrained, Tr7 =Tr7„(o), we have Tr(p^'3' - P^(o) ) log (o) = -3 log(wA^^B) sinh2r(TrY^-(J^V)^ - hn- 
). For Vj = (7j — l)/(7j + 1) < l,thus— \og{v aVb) > 0. What left to be proved is that at the con- 
straint of Tr7 = AnE = 4En, —{Jn^^Y ~ reaches its maximum when 7 = 7„(o) = 2El2n- We start with 

any given 7 with Tr7 = 4i?„. 7 can be symplectically diagonahzcd to SjS'^ = (i?'ag{7i, 71, 72, 72, • • • .7n!7n}- 
The symplectical transformation S can be written as i?2.Di?i, where Ri, R2 are rotations and D = diag{di, l/di, 
c^2,l/c^2 ) dn,l/dn} is the squeezing operation. The rotation does not change the trace of the CM, the 
squeezing operation reduces the trace of the CM in the diagonalizing process. Hence Tr7 > 2^-7^. Let k = 

(Tr7)/(2 J2i tO > 1, 7' = k • rfm5{7i, 71, 72, 72, • • • , 7n, 7n} so that Tr7 = Trj'. We have Tr^ -{Jn^j)^ - hn = 

2J2i y/ii - 1 < Try^ -{Jn^l'Y - hn- For the diagonal CM 7' with the energy constraint Tr7' = 4£„, it is 
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easy to elucidate that when all the diagonal elements are equal, Try — (J„ ^7)^ — = 2 J^i ~ 1 reaches 

its maximal value AE^ — 1. We have proved that 

Trip^.""' -P%nogp"J^<0, (14) 

together with Tr((/9^)— logp'^^p^ = 0. Hence Pn{o) is the extremal state that maximizes CI as far as CI is 
convex at the vicinity of Pn(o)- Similarly, if CI is concave at the at the vicinity of Pn{o), Pn(o) is the extremal state 
that minimizes CI. CI is the difference of two convex function, quantum mutual information and the entropy of 
the input state. This will provide the other way of obtaining the extremal state. 

The next part of this paper is to give evidence of CI really reaching its maximal at Pn{o) if the input energy 
is strong enough. The calculation is based on perturbation theory. Suppose the input state p„ have a complex 
characteristic function Xn{p) =TT{pnD{p)) = Xn(o)(M)(l_+ M*)), where Xn{o){p) = exp[-(7V5 + i) |^|^] is 
the complex characteristic function of Pn(o) i with Ng = E — ^ being the average photon number of the thermal 
input Pi(o) and p = {pi,p2, ■ ■ ■ ,Pn)- The perturbation item can be expanded with power of p and 

p* . The condition that the first and second moments of pn are equal to that of Pn{o) leads to all the item of 
the linear and square power in the expansion of f{p,p*) being 0. The cubic item p^~'^p*'^ as well as other odd 
power items will have no contribution in the first order perturbation. So the first none zero contribution will be 
the fourth power. The other even power items can be neglected comparing with the fourth power as the input 
energy become strong enough. So we suppose 

f{l^,y*) = ^c^j\p^pj\^ . (15) 

It should be noted that items such as P1P2 \Pi\^ and p\p'^ will also contribute to the first order perturbation, 
but they can also be neglected at the strong input energy as we will see below. All other items with unequal 
number of p and p* will not contribute to the first order perturbation. 

Lets first consider the \p\\^ item, the input state now is a direct product of perturbed first mode and 
other n — 1 thermal state modes. The situation is reduced to deal with the perturbation problem of Xi(Mi) — 
Xi(o)(a*i)(1 We have pi = (1 + e J^)Pi(o) = Pi(o) Pi(o) = (1 ~ "s) E^o 1^) (^1 ' ^i^^i 

Vs = Ns/{Ns + 1). The strict eigenvalues of pi are Afc — Xk(a) + £0fc, with Xk{o) = (1 ~ ^s)^^ and ipk — 
Afc(o)[2 — 4:k/Ns + k{k ~ 1)/A^^]. The entropy of p\ can be expanded up to the second derivative as S(^p\) = 
S'(/9i(o)) — 5^^ X]fc '/'fc/-^fe(o) + o(£^), where null first and second moments of are used. The calculation of 
the entropy of p\ is straightforward, it is S{p'y) = S'(p^Q^) — h^'^\- N'{N'+i) ^'^ + o(e^) , with N' = Ng + N. 
The purification of pi(o) is pf^ = J2km V^k^m \kk) {mm\ , such a purification is more frequently used in 
literature but different from our above purification. The state p^ then is expanded in e to the linear item 15:. 

Pi'^ = Pm + ' with $ = i(<i>o + $J), $0 = (1 - ^^.)'[l - 4al6l(l - v,)v7'^^ + a^^il - v,)^vj^]pfgy It 
can be proved that 

{£ <E> I)at/6tVfS = v;^^'a\'a{pfg'^ . (16) 

Thus p"'^' = pf(^)' + £$', with = i($o + $;/), = (1 - Ws)2[2 - 4:a\ai/N, + a\^al/N^]pfl^y Note that 
the trace and all first and second moment of $' are null. The eigenstates of p^^j are Vi \km) with eigenvalues 
Afc„i(o) = (1 — va)va (1 ~ ''^b)vb, where Vi diagonalizes pf^-j and 

ViUiVi = oi coshr — fe| sinhr, (17) 
VibiVi = &i coshr — a| sinhr, 

with tanh2r = 2^/n7{n7TT}/{Ns + N' + 1). The first order perturbation to the eigenvalue will be ^'^m — 
{km\ Vi^'Vi \km) which is 

= Afo„(o)(l-Ws)'{2-4[fccosh2r + (TO + l)sinhV]/iV, (18) 
+ [k{k - l)cosh'*r + (m + l)(m + 2) sinh"^ r + 4fc(m + 1) sinh^ r cosh^ r] /A^f } . 
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Up to item, the entropy wiU be S{pf^ ) ~ S{pf^-^ ) — ^e^ J2km ^km/^km(o) ■ After the summation and taking 
the hmitation of Ng oo, the total increase of CI between the input pn and Pn{o) will be 

^hm^[/,(p„) - /e(p„(o))] = -^^'[^ - 2!^] < 0- 

the positive part which comes from p^'^ state is only | of the negative part which comes from p„. Thus 
Ic{Pn(o)) is maximal in this situation. 

The next perturbation item is \p.ip2^ ■ We only need to deal with the first and second modes with 
X2(pi,P2) = X2(o)(mi7M2)(1 + £ 1/^1/^21^) [T?)j while other modes are kept in thermal states and irrelative. Here 
the degenerate perturbation is applied. The calculation of the entropy difference of P2 and P2(q) is easy because 
the perturbation operator is diagonal in degenerate subspace. The result is 

g(p;)-g(p^(o)) = - 2^,.(^'^,). - (20) 

The calculation of the entropy difference of p^'^ and p^(^| will encounter with non diagonal operators aiflj&i&J 

and a\a2b\b2 in degenerate subspace which indicate inter- mode particle transfer, but in the degenerate subspace 
the total particle number of Q system (or R system) is conserved . The entropy difference can be calculated 
by first summing up in the degenerate subspace then the total particle number of Q system and R system. In 
the summation TrM^^^^^ is involved, where M' is the perturbation operator in the Ith degenerate subspace, 
fortunately it is X^feimi -^feimi special structure of A/'. The final result after the summation and taking 

the limitation of A^s ^ 00 will simply be 

^Iim^^(pr)-^(PS')--Y^- (21) 
Still it is I of the entropy difference of ^^^d Pj^q)- Thus we have 

^lim^[^c(Pn) - /c(p„(o))] = < 0. (22) 

In the situation of X2(pi, P2) = X2(o) (mi; M2)[1 + ^(Ipi |^ + c |pip2|^)] , the entropy difference will be the sum 
of each term because the cross item is null by the null of the first and second moment of and <&'. The general 
case of /(p, p*) will be 

^nm^[/,(p„) - /,(p„(o))] = -^^(4^4 + E 4) < 0- (23) 
The conclusion is that thermal state of infinitive energy achieves the maximal of coherent information. That 

is 

lim maxi/c(p„,£®") = max{0,-log2(e7V)}, (24) 

{e — 2.71828 . . .).The result is verified up to the nonzero lowest power of the inverse of the input energy in each 
channel use. Whether it is correct for the state without an item contributing to N~'^ is not known. 

We obtain a local maximum, whether it is the global supremum should be verified. For lower noise channel, 
this can be verified. For a given channel, the CI difference is a function of Ng. In the two case we calculated, 
the CI difference is positive infinitive at iV^ 0, as Ng increases, it monotonically decreases to 0. After that it 
decreases further to negative then increases but still keeps negative and never turns to positive. At Na 00, it is 
negative as we elucidated above. Denote the zero point as N^o, then calculate the quantum mutual information 
I{Pn{o){Nso)i N) which is the supremum of all state with equal or less energy. I{pn,N) is greater Ic{pn,N). 
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Hence if /(pi(o) (iVso), ^) < ~log2(e^) : we have CI to be less than — rilog2(eiV) in the interval [0, iVso], thus 
for N < Nc we can safely conclude that 

Q = max{0,-log2(e7V)}, (25) 

as far as the high order perturbation and high power of do not destroy the maximal property of CI at infinitive 
input energy, where Nc is the solution of /(pi(o) (iVso), Nc) = — log2(eiVc) . We have Nc = 0.1756 which comes 
from the two-mode perturbation. 

For the n uses of the thermal noise channel with TV < N^ the supremum of the whole CI is achieve by an 
input of the direct product of the identical thermal noise states. It is followed that the quantum capacity of 
thermal noise channel is equal to the one-shot quantum capacity dlJI] of the channel. The achievable of the 
quantum capacity of the thermal noise channel by quantum error-correction codes had been proven 1 101. 
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